PERFECT FLUID DARK MATTER 



F. RahamarQ 

Department of Mathematics, Jadavpur University, Kolkata 700 032, West Bengal, India 

K. K. Nandi 

Department of Mathematics, North Bengal University, Siliguri - 734013, West Bengal, India 

A. Bhadra^ 

High Energy and Cosmic Ray Research Centre, University of North Bengal, Siliguri, WB 734013, India 

M. Kalani 

Department of Physics, Netaji Nagar College for Women, Regent Estate, Kolkata-7 00092, India 

K. Chakrabort>§ 

Department of Physics, Government Training College, Hooghly - 712103, West Bengal, India 

(Dated: September 23, 2010) 

Taking the flat rotation curve as input and treating the matter content in the galactic halo region as perfect fluid, 
we obtain space time metric at the galactic halo region in the framework of general relativity. We find that the 
resultant metric is a non-relativistic-dark-matter-induced space-time embedded in a static Friedmann-Lemaitre- 
Robertson-Walker universe. This means that the flat rotation curve not only leads to the existence of dark 
matter but also suggests the background geometry of the universe. The flat rotation curve and the demand that 
the dark matter be non-exotic together indicate a (nearly) flat universe as favored by the modern cosmological 
observations. We obtain the expressions for energy density and pressure of dark matter halo and consequently 
the equation of state of dark matter. Various other aspects of the solution are also analyzed. 

PACS numbers: 04.20 Gz,04.50 + h, 04.20 Jb 
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I. INTRODUCTION 

It has been known for a long time that rotation curves of 
neutral hydrogen clouds in the outer regions of galaxies can- 
not be explained in terms of the luminous matter content of 
the galaxies, at least within the context of Newtonian gravity 
or of general relativity (see [1] for a review). In addition, the 
velocity dispersion in galaxy clusters indicates a much greater 
mass than would be inferred from luminous matter contained 
in the individual galaxies. The most widely accepted expla- 
nation, based on the standard gravitational theory, postulates 
that almost every galaxy hosts a large amount of nonluminous 
matter, the so called gravitational dark matter, consisting of 
unknown particles not included in the particle standard model, 
forming a halo around the galaxy. The dark matter provides 
the needed gravitational field and the mass required to match 
the observed galactic flat rotation curves in galaxy clusters. 

The dark matter problem arises because of a naive New- 
tonian analysis of the near constant tangential velocity of ro- 
tation up to distances far beyond the luminous radius of the 
galaxies. This leads to the conclusion that the energy density 
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decreases with the distance as r 2 and therefore that the mass 
of galaxies increases as m(r) oc r. 

There are several proposals for the dark matter component, 
ranging from new exotic particles such as those predicted by 
supersymmetry [2] to other less exotic candidates such as 
massive neutrinos or even ordinary celestial bodies such as 
Jupiter like objects. Several other analytic halo models ex- 
ist in the literature including those sourced by scalar fields. 
Fay [3] considered modelling by Brans-Dicke massless scalar 
field while Matos, Guzman and Nunez [4] considered mass- 
less minimally coupled scalar field with a potential. A bo- 
son star formed by a self interacting massive scalar field with 
quartic interaction potential as a model of galactic halo was 
investigated by Colpi, Shapiro and Wasserman [5]. A similar 
boson star as a model of galactic halo was first investigated 
by Lee and Koh [6]. A recent halo model is discussed also 
in the braneworld theory [7]. Some authors [8] have consid- 
ered global monopoles as a candidate for galactic halo in the 
framework of the scalar tensor theory of gravity. 

Alternatively, explanation of the observed flat rotation 
curve is provided through modification of Newtonian dy- 
namics (in the region of very small accelerations) [9,10]. 
Another attempt to resolve the dark matter problem lies 
within the framework of Weyl conformal gravity suggested 
by Mannheim and Kazanas [11], with the distinction that it 
does not postulate flat rotation curve as an input but predict it. 
Such a prediction is possible also within the Chern Simmons 
gravity inspired by string theory [12,13]. 
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In recent years, the amount of dark matter in the Universe 
has become known more precisely: CMB anisotropy data in- 
dicate that about 85% of total matter in a galaxy is dark in 
nature. Big bang nucleosynthesis and some other cosmolog- 
ical observations require that the bulk of dark matter be non- 
baryonic, cold or warm, stable or long-lived and not interact- 
ing with visible matter. However, despite long and intensive 
investigations, little is as yet known about the nature of dark 
matter. 

The purpose of the present work is to show that, with the 
input of flat rotation curve and the assumption that dark mat- 
ter be described as a perfect fluid, the general theory of rel- 
ativity not only accounts for the dark matter component of 
galactic clusters but also suggests the spatial curvature of the 
universe. The gravitational influence of an arbitrary dark mat- 
ter component is controlled by its stress tensor. In this con- 
text we note that, while the existing approaches of explain- 
ing flat rotation curves are useful in their own right, many 
of the models end up with predicting anisotropic dark matter 
fluid stress tensor. On the other hand, no physical mechanism 
is stated explaining why a spherical distribution should have 
such anisotropy. There is also no observational evidence in 
support of it. Therefore, it seems more reasonable to consider 
an isotropic perfect fluid distribution for dark matter because 
predictions from such model at stellar and cosmic scales have 
been observationally corroborated beyond any doubt. 

We shall particularly study the general features of dark mat- 
ter such as its equation of state. The mass distribution in galac- 
tic haloes should provide a direct probe into the nature of par- 
ticles constituting the dark matter because the inner structure 
of the halo is particularly sensitive to the dark matter proper- 
ties [14]. In this investigation, we shall work only from a fluid 
perspective leaving the question about the particle identity of 
dark matter open. 



II. GRAVITATIONAL FIELD IN THE DARK MATTER 
REGION 



where the functions v(r) and A(r) are the metric potentials. 
Then the Einstein field equations become (c = 1): 
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For a circular stable geodesic motion in the equatorial plane, 
the consideration of flat rotation curve gives the condition (see 
Appendix) 
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where I is given by I — 2(w^) 2 and B is an integration con- 
stant. The observed rotational curve profile in the region dom- 
inated by dark matter is such that the rotational velocity be- 
comes approximately a constant with <~ 10~ 3 ( 300km/s) 
for a typical galaxy. The Eqs. (3) and (4) then lead to the 
following equation 
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Galactic rotational velocity profiles [15] of almost all the 
spiral galaxies are characterized by a rapid increase from the 
galactic center, reaching a nearly constant velocity from the 
nearby region of the galaxy far out to the halo region. Our 
target is to exploit this observed feature to obtain the space- 
time metric in the halo region, and analyze it. As mentioned 
already, observational data suggest that the dark matter com- 
ponent in the galaxy accounts for almost 85% of its total 
mass. Naturally, luminous matter does not contribute signif- 
icantly to the total energy density of the galaxy, particularly 
in the halo region. Therefore, we shall treat the matter in 
the galactic halo region as a perfect fluid defined by stresses 
77 = Tg = = p, where is the matter energy momen- 
tum tensor. 

The general static spherically symmetric spacetime is rep- 
resented by the following metric 



ds 2 = -e» (r Ut 2 + e x{r Ur 2 + r 2 {dd 2 + sin 2 9d^ 2 ), (1) 
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An exact solution of Eq.(6) is given by 
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where D is an integration constant. Several observations fol- 
low from the metric (1): 

(a) Note that the space time metric given by Eq.(l) through 
Eqs. (5) and (9) is an interior solution. This type of spacetime 
definitely cannot be asymptotically flat neither can it have the 
form of a spacetime due to a centrally symmetric black hole. 
What can be said is that this line element describes the region 
where the tangential velocity of the test particles is constant 
and that it has to be joined with the exterior region with other 
types of space time. The extent of v^= constant ends at some 
larger distance where the region becomes asymptotically flat. 
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In principle, the constant D should be obtained from the junc- 
tion conditions but the galactic boundary is not observation- 
ally defined yet. It is more likely that dark matter distribution 
will not remain the same near the periphery of the halo. 

(b) Inserting the above solution in Eqs. (2) - (4), one can 
readily get the expressions for p, p as 
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The most notable aspect is the presence of the last term 
in the expressions of both energy density p and pressure p 
that increases with radial distance (from galactic center). Note 
that according to the Newtonian theory, which is supposed 
to be indistinguishable from general relativity in very weak 
field, one should expect only the term pNewton = sttGj^ m 
agreement with the Poisson equation, which integrates to give 
the Newtonian mass M(r) increasing linearly with r. But 
the input of the constant tangential velocity leads to general 
relativistic corrections to Newtonian expressions as evident in 
Eqs. (10), (11). The first term of the right hand side of Eq.( 10) 
gives the expected Newtonian term in the leading order and to 
the same order a general relativistic correction term ^V^-Pr. 
Since I is small, the correction terms are small as expected. 

However, the second term in Eq.(10) has completely differ- 
ent nature than the conventional feature of dark matter energy 
density; the corresponding mass increases nonlinearly with 
radial distance r, though a bit slowly as I is very small. This 
contribution will vanish if D vanishes. We would see that even 
D = is consistent with the flat rotation curve. 

(c) The parameter D is recognized [17] as the spatial cur- 
vature (with a negative sign) of the universe. A comparison 
of the obtained space time metric in the limit I -> with the 
static Friedmann-Lemaitre-Robertson- Walker (FLRW) met- 
ric [ ds 2 = -dt 2 + Y^idr 2 + r 2 (d0 2 + sin 2 6d<j) 2 ) ] leads 
to the identification that D = —k. Since our input was only 
the flat rotation curve and we have not considered anything 
regarding cosmological spatial curvature while deriving the 
metric, the appearance of D in the metric element is quite in- 
teresting. The solution thus may be thought of dark matter 
induced space-time embedded in a static FLRW metric. In 
general I ^ and hence a is not equal to —2, and the inter- 
pretation of spatial curvature does not seem at all evident. 

An immediate question is how a spatially curved FLRW 
universe could be static. If it is due to balance of the curvature 
term by some other fluid with equation of state p = —3p, the 
question of stability will arise. It seems that we have ended 
up with a static universe because right from the beginning we 
looked for a static solution (all the metric element are consid- 
ered time independent). While working on a local problem 
(flat rotation curve), the universe is usually considered at any 
particular epoch fixing the scale factor R(t) to be constant (of- 
ten normalized to unity, R(t ) — 1 at the present epoch). The 
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FIG. 1 : The variation of uj with D for a typical distance r — 200 in 
Kpc. 



dynamicity of the FLRW metric is thus absent in local gravita- 
tional phenomena. However, we see in the present work that 
the curvature effect has appeared even in local gravitational 
phenomena. It would certainly be an interesting but challeng- 
ing task to study the stability criterion in this configuration 
[17]. We reserve it for a separate communication. 

(d) The equation of state parameter cj for the effective fluid 
(dark matter plus 'curvature fluid '), can be obtained directly 
from Eqs. (10) and (11), which is given by 
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Since I is small, effectively u> ~ j^rrETD ■ ^ tne tota l mat " 
ter in the flat rotation curve region of galaxies has to be non- 
exotic, i.e. if the dark matter satisfies the known energy con- 
ditions, oj for the total matter content must be positive. This 
is achieved only if we set l 2 r a /4 < D < lr a /3. Note that 
our discussion is restricted only in the region of flat rotation 
curve of galaxies for which r is typically between few ten kpc 
to few hundred kpc and as mentioned already is few hun- 
dred km/s. Hence positive uj implies that D is nearly zero, if 
not exactly. Figure 1 shows that variation of u with D for a 
typical distance r = 200 kpc and — 300 km/s. 

As may be seen from the figure 1, the limit on D is much 
stringent compared to the (cosmological) observational re- 
striction. The equation of state of the dark matter component 
may be obtained by taking D = in the above equation which 
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gives 
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Thus p << p for dark matter implying its non-relativistic na- 
ture, which is a well known fact. Further u is positive, which 
means that the fluid is non-exotic. On the other hand, the 
equation of state of 'curvature fluid 'is found to be p + 3p = 
from Eq.(12) in the limit I — > 0, which is a familiar result. 
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which is often convenient. Such a form has the advantage that 
it immediately reveals not only the mass parameter m(r) but 
also shows that the proper radial length is larger than the Eu- 
clidean length because r > 2m(r). But most importantly, this 
inequality is essential for signature protection, which dictates 
that e A > 1. This is a crucial condition to be satisfied by any 
valid metric. Now, from the metric function (9), for D = 0, 
we get 



e A = H : > 1. 
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This implies that the essential requirement is fulfilled for this 
value of D. 

(f) Let us rewrite the metric (1) for D = under the radial 
rescaling 



which yields 
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FIG. 2: The variation of R with r in Kpc for small different values 
of D. We choose ~ 10 -3 (300km/s) for a typical galaxy. 



grows as ~ 10 in units of flat surface area, which indicates 
another deviation from the massless scalar model in which it 
grows as ~ 10~ 12 [4]. 

(g) The Ricci scalar for the derived spacetime is given by 
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As I — > 0, R = —6D, once again suggesting that D is the 
spatial curvature. We plot R vs r. One may note that the value 
of R is small. For small different values of D, one can not 
distinguish the variations of R with respect to r ( see figure 2 
). But for higher values of D, figure 3 shows a sharp variation 
of R. Note that large values of D imply cj(r) < 0, i.e., dark 
matter has to be exotic in nature. 
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showing a surplus angle in the surface area given by 
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HI. OTHER ASPECTS 



The total gravitational energy: 
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If the probe particles were photons, so that = 1, the surface 
area would remain finite but reduced to half of the spherical 
surface area S2 — Airr 2 . This is an interesting result, which 
distinguishes itself from that in the massless scalar field model 
where S — > 00 as vr = 1 [4]. For a typical rotational velocity 
of = 10 -3 in the galactic halo region, the difference of the 
two surface areas 
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One notes from equation (12) that the halo matter is not 
exotic in nature and consequently, we expect attractive grav- 
ity in the halo. Following the suggestion given by Lyndell - 
Bell et al [18] , we calculate the total gravitational energy Eq 
between two fixed radii , say, r-y and r2'- 
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FIG. 3: The variation of 7? with r in Kpc for large different values 
of D. We choose ~ 10~ 3 (300km/s) for a typical galaxy. 
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is the Newtonian mass given by 



M — 47r / pr dr 



(23) 



FIG. 4: The variation of Eg with r in Kpc. The lower limit of 
integration in equation (24) fixed at, say, n — 100 kpc while V2 is 
varied from 100 to 500 Kpc. We choose G = 1, D = 10 -11 and 
v <t> ~ 10~ 3 ( 300km/s) for a typical galaxy. 
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Thus we get the total gravitational energy as 
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The figures 4 and 5 show that the total gravitational energy 
is small but negative whether we choose D non-zero or zero 
for arbitrary T2 > r\ > 0. Thus for the existence of non- 
exotic matter in the halo, we recommend the value of D < 
10 -11 . In the distant halo region we have taken, typically , r 
~ 200 kpc in the figures below. 
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Attraction: 



Now we study the geodesic equation given by 



dV* dx»dx< 
dr 2 a dT dr 
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FIG. 5: The variation of Eg with r in Kpc. The lower limit of 
integration in equation (24) fixed at, say, ri = 100 kpc while r2 
is varied from 100 to 500 Kpc. We choose G = 1 , D = and 
v <p ~ 10~ 3 ( 300km/s) for a typical galaxy. 
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for a test particle that has been "placed' at some radius r . 
This yields the radial equation 
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which is negative as the the quantity in the square bracket is 
positive. Thus particles are attracted towards the center. This 
result is in agreement with the observations i.e., gravity on the 
galactic scale is attractive ( clustering , structure formation etc 

)• 



C. Stability: 



Let us define the four velocity U a = for a test particle 
moving solely in the space of the halo (restricting ourselves 
to 9 = it/2), the equation g ucr U' / U' J = —to 2 , can be cast in a 
Newtonian form 



which gives 
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IV. CONCLUSIONS: 

It is well known that observation of flat rotation curve sug- 
gests that a substantial amount of non-luminous dark matter 
is hidden in the galactic halo. Here we have found that the 
flat rotation curve suggests also the background geometry of 
the universe. The space-time geometry we have obtained can 
be interpreted as the one due to dark matter embedded in the 
=r ° static FLRW universe. This is probably the first indication 
that the spatial curvature of the universe can be obtained from 
a local gravitational phenomenon. 

If we demand that matter in the flat rotation curve region 
be non-exotic (i.e., obey the usual energy conditions), we ob- 
tain the result that the universe should be nearly flat, if not 
exactly so, which is consistent with modern cosmological ob- 
servations. 

The equation of state of the dark matter component has 
been obtained by treating it as perfect fluid and the expres- 
sions for general relativistic correction terms for the pressure 
and energy density over those obtained from the Newtonian 
theory are also derived. The corrections are, however, small 
as expected. 

The geodesic equation [Eq.(26)] of a test particle in the de- 
rived spacetime suggests that the particle will be attracted to- 
wards the center. We have quantified the attractive effect in 
the relativistic case by calculating the total gravitational en- 
ergy Eg ( which is negative ) in the halo region. We have also 
demonstrated the stability of the circular orbits in our space- 
time solution. Thus our solution satisfies two crucial physical 
requirements - stability of circular orbits and attractive gravity 
in the halo region. Investigation on other observational con- 
straints on the model is underway. 
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where the constants E and L, respectively, are the conserved 
relativistic energy and angular momentum per unit rest mass 
of the test particle. Circular orbits are defined by r — 
R =constant. so that §f = and, additionally, \ r =R= 0. 
From these two conditions follow the conserved parameters: 
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and using it in V(R) = —E 2 , we get 
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The orbits will be stable if \ r =R< and unstable if 
-r> 0. Putting the expressions for L and E in 



dr 2 

~j~t \r=R, we obtain, after straightforward calculations, the 
final result, viz., 
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Appendix 

To derive tangential velocity of circular orbits, we start with 
the line element 

ds 2 = -e» (r Ut 2 + e x(r Ur 2 + r 2 (d9 2 + sin 2 6d<i> 2 ). 

The Lagrangian for a test particle reads 

2L = -e^ r H 2 + e A «r 2 + r 2 (6 2 + sm 2 ^ 2 ). 



(33) 

One may note that, for D = as well as for D ^ 0, 
4pr \r=R< 0. So the circular orbits are always stable. 



We guess the conserved quantities, the energy E = e v ^ r H, the 
(^-momentum = r 2 sin 2 9<f>, and the total angular momen- 
tum, L 2 = L e 2 + (Jj£^ , with L g = r 2 6. The radial motion 
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equation can be written as: 

r 2 + V(r) = 
with the potential V(r) given y 

V(r) = -c-W (e"^ 2 " ^ " l) ■ 

For circular orbits, we have the conditions, r = 0, V r = 
and > 0. These imply the following expressions for the 
energy and total momentum of the particles in circular orbits: 
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The second derivative of the potential evaluated at the extrema 



V(r) r 
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Now, the tangential velocity 

(^) 2 = r 2 e -'^)(0 2 + S m 2 ^ 2 ) 
can be obtained for particles in stable circular orbits as 



r(e^ r >) r 
2e"( r ) 



Using the tangential velocity to be constant for several radii, 
the above expression yields 

e v = B a r l 

where I is given by I = 2(v^) 2 and B is an integration con- 
stant. 
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